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Using Monte Carlo simulations, we explore the magnetic phase diagram of the triangular spin 
tubes coupled with a ferromagnetic inter-tube interaction for CsCrF 4 . A planar structure of the 
coupled tubes is topologically equivalent to the Kagome-triangular lattice, which induces nontrivial 
frustration effects in the system. We particularly find that, depending on the inter-tube coupling, 
various ordered phases are actually realized, such as incommensurate order, ferromagnetic order, 
and Cuboc order, which is characterized by the non-coplanar spin structure of the twelve sublattice 
accompanying the spin chirality breaking. We also discuss a relevance of the results to recent 
experiments of CsCrF 4 . 
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I. INTRODUCTION 

Recently, a triangular spin tube has attracted much 
interest, where its geometrical frustration and quasi-one- 
dimensionality cooperatively induce exotic magnetic be¬ 
haviors. Indeed, theoretical investigations of the ^ = 1/2 
quantum spin tubei have revealed various interesting 
properties such as gapful ground state^r— , field induced 
chirality orderiSiiS, etc. Moreover, extensive theoreti¬ 
cal researches have been performed for various quantum 
spin tubes, such as integer-spin tubesii, S = 3/2 tri¬ 
angular tubesi^d^, and four-leg tube^d. Also, the tri¬ 
angular spin tube has been a target of intensive exper¬ 
imental studies. For example, several experiments on 
[(CuCl 2 tachH) 3 Cl]Cl 2 , which is a S' = 1/2 spin tube 
consisting of alternating triangles along the tube direc¬ 
tion, clarified various characteristic behaviors originat¬ 
ing from the tube structur o^^d^ . Moreover, straight-type 
spin tubes CsCrF 4 and a-KCrF 4 have been recently syn¬ 
thesized, which are respectively based on equilateral and 
non-equilateral triangles^^^— . These compounds inter¬ 
estingly provide essential information about the shape 
dependence of the unit triangle in the spin tubes. 

In CsCrF 4 , Cr^+ ions having S = 3/2 spin form a rigid 
equilateral triangular tube (Fig. 1), where dominant ex¬ 
change couplings are antiferromagnetic and an inter-tube 
coupling is estimated to be basically very small. Accord¬ 
ingly, no anomaly associated with a phase transition was 
observed by bulk measurements and ESR experiments 
down to T = 1.5Kpd^“— On the other hand, a recent ex¬ 
periment of AC susceptibility observed anomalous slow 
dynamics suggesting a magnetic long-range order below 
4K2i. In addition, a very recent neutron diffraction ex¬ 
periment suggests that this magnetic order is inconsistent 
with a naive 120° structure due to the conventional tri¬ 
angle latticei^ Thus, it is expected that the equilateral- 
triangle structure and a small but finite inter-tube cou¬ 
pling cooperatively induce a non-trivial magnetic struc¬ 
ture in CsCrF 4 , which could be indeterminate in the bulk 


quantities. 

In order to analyze the magnetic structure of CsCrF 4 , 
a key observation is that Cr^+ has a relatively large 
spin S = 3/2, and a certain spin order is suggested by 
magnetic diffraction peaks in the neutron experimenliS^. 
Thus, we can expect that the magnetic order of CsCrF 4 
is basically described by the classical Heisenberg model 
defined on the triangular tube lattice. As will be de¬ 
picted in Figs. 1 and 2, moreover, the inter-tube cou¬ 
pling in the afc-plane has the same lattice topology as the 
Kagome-triangular lattice^ although the exchange cou¬ 
pling along the c-axis is dominant in the spin tube. Thus, 
the lattice structure of CsCrF 4 involves the frustration 
effect even for the ferromagnetic inter-tube coupling. 

Of course, the Heisenberg model on the planar lattice 
has no magnetic long-range order at a finite tempera¬ 
ture. In the present coupled tubes, however, the three- 
dimensional(3D) couplings possibly stabilize a peculiar 
spin fluctuation originating from the Kagome-triangular 
structure. Using Monte Carlo (MC) simulations, in this 
paper, we investigate finite-temperature phase transi¬ 
tions of the classical Heisenberg model on the triangu¬ 
lar spin tubes with the inter-tube interaction. In par¬ 
ticular, we find that the twelve-sublattice spin structure 
—Cuboc order— in the a6-plane emerges in the small 
ferromagnetic inter-tube coupling regime. The Cuboc or¬ 
der, which was originally proposed for the ground state of 
the planar Kagome lattice model with the next-nearest- 
neighbor interaction, is characterized by a non-coplanar 
spin structure with the triple-q wave vectors^i^— . In the 
spin tubes, this Cuboc order can be stabilized to be the 
3D long-range order by the strong leg coupling of the 
tube at a finite temperature. We also discuss nature of 
the transitions for the Cuboc phase, as well as incom¬ 
mensurate and ferromagnetic phases, depending on the 
inter-tube coupling. Finally we discuss the relevance to 
the CsCrF 4 experiments. 

This paper is organized as follows. In Sec. H, we ex¬ 
plain details of the model and the possible orders. In Sec. 
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FIG. 1 : (Color online) Lattice structure of triangular tubes 
with an inter-tube interaction, (a) A triangular spin tube, 
where Ji is the dominant antiferromagnetic coupling in the 
tube direction (c-axis) and J2 denotes the antiferromagnetic 
interaction in the unit triangle, (b) The lattice structure in 
the afe-plane. The triangles of the solid lines correspond to the 
J2 coupling in the spin tubes and the triangles of the broken 
lines denote the ferromagnetic inter-tube interaction J3. 


J2 



FIG. 2 : (Color online) (a) Topology of the coupled spin tube 
lattice in the afo-plane is equivalent to the Kagome-triangular 
lattice, (b) The ground state phase diagram of the classi¬ 
cal Heisenberg model on the coupled triangular spin tubes. 
The left-going arrow with a broken line shows the direction of 
J3 parameter corresponding to the finite-temperature phase 
diagram in Sec. IV. 


Ill, we describe details of MC simulations and definitions 
of order parameters. In Sec. IV, we present results of MC 
simulations and summarize the phase diagram with re¬ 
spect to the inter-tube coupling. We also mention univer¬ 
salities of the phase transitions. In Sec. V, we summarize 
the conclusion and discuss the relevance to the CsCrF4 
experiments. 

II. MODEL AND ORDERS 

As in depicted in Fig. [I] a bundle of triangular spin 
tubes in CsCrF4 runs in the c-axis direction and these 
tubes with the inter-tube coupling cover the triangular 
lattice in the a6-plane. We thus consider the classical 
Heisenberg model on the stacked triangular lattice, which 
reads 

n = JiY. S{r,)-S{rk) + J2 Y. 

+J3 ^ s(n) . S{r,) ( 1 ) 

{(id)) 

where represents the position vector of site i, and 
S{ri) = Sfe^ + S^By + with \Si\ = 1 denotes the 
vector spin at the *th site. Note that {a € x, y, z) 
indicates the unit vector in the spin space, while the 
primitive lattice translation vectors are represented as 
a,b and c with |a| = |b| = |c| = 1 (Fig. [ 1 }. Moreover, 
{i,k)c denotes sum for the nearest-neighbor spins along 
the c-axis, (i, j)a indicates sum of spin pairs in the unit 
triangle, and runs over spin pairs of the inter-tube 

couplings in the a6-plane. In this paper, we basically as¬ 
sume the antiferromagnetic interactions Ji > J2 > 0 in 
the spin tube and the ferromagnetic inter-tube interac¬ 
tion J3 < 0 . Note that the LDA-I-U calculation gives 
J1/J2 ~ 2.0 and J3/J2 — 0 with J2 ~ 20K for CsCrF^Sl. 

In analyzing possible ordering of the coupled-tube 
model, an important point is that the dominant cou¬ 


pling Ji along the tube direction does not cause any 
frustration. Thus, the staggered pattern of the spin order 
formed in the a6-plane is realized in the c-axis direction. 
This implies that the low-temperature spin structure is 
essentially attributed to the frustrating interactions in 
the a6-plane, which we will actually justify with MC sim¬ 
ulations in the next section. In the following, we therefore 
assume the staggered order in the c-axis direction, and 
concentrate on the spin structure in the oti-plane. 

As in Fig. UKa), the planar structure of the model 
is topologically equivalent to the Kagome-triangular lat¬ 
tice, where J3 corresponds to the nearest-neighbor cou¬ 
pling competing with the next-nearest-neighbor interac¬ 
tion Then, a candidate of the ground-state order is 
classified by Fourier transformation of the Hamiltonian. 
Defining the unit cell as a triangle of the J2 coupling, we 
have 

^ = 2 Yl ( 2 ) 

where Sq^a = Sr Here, r represents 

the position of a unit triangle, a G {A, H, C} indicates 
the sublattice index in the unit triangle, and N is the 
total number of the unit triangles in the system. In ad¬ 
dition, Japiq) = Sa,/3 where r„/3 and Jajs 

respectively denote the relative vector of a spin pair and 
the corresponding coupling associated with spins in the 
unit cell. Note that the wavevector q runs over the 1 st 
Brillouin zone in the ab plane. By determining the lowest 
energy state of Jap{q), we have the ground-state phase 
diagram of Eq. (P) in the J2-J3 plane (FigPb)), which 
is equivalent to that obtained in Ref. [ 23 j |. 

In Fig. HKb), CsCrF4 is located nearby J3 ~ 0 and 
J2 > 0 . If J3 is antiferromagnetic, Eq. © is a triangular 
lattice antiferromagnet for which the ground state is the 
120 ° structure. However, the neutron diffraction experi¬ 
ment suggested that the order of CsCrF4 is not explained 
by a naive 120 ° structure^^. We thus discuss the negative 
J3(< 0) region, where the nontrivial exotic order phase 































3 


actually emerges; In — J3 < J2 region, particularly, the 
minimum of Japisi) located at M point in the Brillouin 
zone, where the non-coplanar order with the twelve sub¬ 
lattice that is called ’’Cuboc” order can be stabilized. As 
—J3 increases, the incommensurate order appears in the 
J2 < —J3 < 2J2 region, and finally the ferromagnetic 
order becomes stable for — J3 > 2J2. This ground-state 
phase diagram suggests that, in CsCrF4, the Cuboc state 
can be stabilized to be the 3 D long-range order at a fi¬ 
nite temperature by the tube-leg coupling Ji, even if the 
amplitude of the inter-tube interaction J3 is very small. 

Here, we briefly summarize the essential property of 
the Cuboc order in 0 < — J3 < J2, which was orig¬ 
inally found in the Ji — J2 Heisenberg model on the 
Kagome lattice^^—. The Cuboc order has a non- 
coplanar spin structure accompanying the spontaneous 
symmetry breaking of the lattice translation. Fig. Ha) 
shows the extended unit cell for the Cuboc order in the 
a&-plane, where we assign a number to each of four tri¬ 
angles, for later convenience. Then, the strong antifer¬ 
romagnetic coupling J2 basically imposes the 120° struc¬ 
ture in each triangle. Then, a significant point is that 
the 120°-structure planes can relatively tilt among the 
four triangles, so as to reduce the energy due to the frus¬ 
trating J3 interaction. Gluing the four tilting triangles 
of the 120° structure, we obtain the tetrahedron where 
the twelve spins in the extended unit cell are attached. 
As shown in Fig. Hb), the three spins on the unit tri¬ 
angle of the original lattice are mapped into the vertices 
of the corresponding triangle on the tetrahedron, where 
the vector-spin chiralities sitting on the four 120° planes 
point the radial direction from center of the tetrahedron. 
In this sense, the vector-spin chirality associated with the 
tetrahedron can be a good order parameter of the non- 
coplanar Cuboc spin structure. Note that, if the spin 
vectors are arranged at the origin of the spi n sp ace, we 
have the same schematic diagram as in Refs. j 24 l - [^ . 

For the Cuboc order, the magnetic propagation vectors 
interestingly have a triple q structure in the aAplane, 
reflecting the above characteristic spin configuration. Let 
us write the sublattice spin in the unit cell at position r as 
Sa{r), where a G {A, B, C} is the sublattice index in the 
unit triangle. Then, the Cuboc order can be explicitly 
written as 

SA(r) = cos - r) - cos ^ 

SB(r)=cos(i,,.r)^-cosQ,..r)^ 

Sc{r) = cos - cos (^Qa ■ ( 3 ) 

where Qa , are the reciprocal lattice vectors of the prim¬ 
itive translation vectors a, b, and = qa ~ qb- This 
triple q is an important feature of the Cuboc phase, 
and plays an essential role in the analysis of the neu¬ 
tron diffraction experiment. Finally, we note that, in the 
c-axis direction, the magnetic propagation vector is <7c/2. 



FIG. 3: (Color online) (a) The extended unit cell for the 
Cuboc order with the twelve sublattice structure. The label 
{A, B, C} indicates each vertex in the unit triangle and the 
number {1,2,3,4} represents the label of a triangle in the ex¬ 
tended unit cell, (b) The arrow with the numbers shows the 
spin directions forming 120° structure in the corresponding 
unit triangle. The four tilting 120°-structures in the extended 
unit cell form the tetrahedron represented by the dashed lines. 


where the simple staggered pattern appears. 


III. MONTE CARLO SIMULATION 

In the previous section, we have discussed the ground- 
state orders of the triangular spin tubes in the aft- 
plane for — T3 < J2. So far, investigations of the pla¬ 
nar Kagome-triangular-lattice model have clarified that, 
although there exists no true long-range order of the 
spin, the chirality degrees of freedom associated with the 
Cuboc order induces the Z2 symmetry breaking transi¬ 
tion even at a finite temperature^ii^^. For the coupled- 
tube system, which contains the full 3 D couplings, we 
can expect the finite temperature transitions associated 
the Cuboc long-range order, the incommensurate order, 
as well as the ferromagnetic order. In order to address 
the nature of these finite-temperature transitions, we per¬ 
form extensive Monte Carlo simulations for the coupled- 
tube system of Eq. © 

Here, we comment on the notation of the system size. 
In the following, we basically represent the linear dimen¬ 
sions of the system by the number of triangles associated 
with J2 couplings. Thus, La{Lb) means a number of tri¬ 
angles in the a(&)-axis direction, and Lc denotes length of 
a tube in the c-axis direction. In this paper, we basically 
deal with the system of La = Lb = Lc = L, for which the 
total number of spins in the system is N = 3 L^. 


A. Details of simulation 

We employ the Wolff’s cluster algorithn^ combined 
with the Metropolis local update. Usually, the Wolff al¬ 
gorithm is not efficient for frustrated systems, where a 
large cluster containing almost all of spins are often gen- 
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erated. However, we find that the coupled spin tubes 
has no frustration in the c-axis, which makes possible 
cluster growing of an efficient size in the tube-direction. 
Thus the Wolff algorithm works very well for the coupled 
spin tubes. Note that the parallel tempering method^^ is 
additionally used in practical computations, if the relax¬ 
ation to the equilibrium is difficult. On the basis of the 
above algorithm, we have performed MC simulations for 
the system of L = 8,16, • • • ,36 with the periodic bound¬ 
ary condition. We particularly explore the T — J 3 phase 
diagram along the broken line in Fig. [2Ib) with J 2 = 1.0 
fixed. Typical numbers of MC samples are 2^® ~ 2^^. 


B. Order parameters 


As pointed out in the ground-state phase diagram of 
Fig. [2jb), the coupled triangular spin tubes have the 
various ordered states. To classify these ordered states 
in simulations, we need to define appropriate order pa¬ 
rameters. In particular, the Cuboc order has the non- 
coplanar spin structure with the twelve sublattice in the 
a 6 -plane, which is staggeredly stacked in the c-axis di¬ 
rection. Taking account of this structure, we define the 
Cuboc sublattice magnetization rUa^p in a certain ab- 
plane, where a £ {A, H, C} indicates a vertex of a trian¬ 
gle, and /3 £ {1,2,3,4} specihes a triangle in the extended 
unit cell (Fig. EJa)). 


TTloi,0 


1 


rp 


( 4 ) 


where rp denotes the position of the triangle labeled by 
/3 in the extended unit cell, and Ylirp represents the sum 
with respect to triangles having the same a in the plane. 
In addition, the normalization = L^/4 is the num¬ 
ber of the extended unit cell in the ab plane. Then, the 
sublattice magnetization order parameter in the entire 
system is defined as 

= (5) 


where ic denotes index of the tube direction (Fig. HJa)). 

In — J 3 > 2 J 2 region, the ground state of the system 
is the ferromagnetic ordered state in a certain Kagome- 
triangular plane . Thus, we also introduce the ferromag¬ 
netic order parameter as 

( 6 ) 

a .,0 


which detects the uniform magnetization in a Kagome- 
triangular plane. Then, this in-plane ferromagnetic order 
is staggeredly stacked in the c-axis direction. Thus we 
define the total ferromagnetic order parameter Mp as 

MF^i^(-ir‘=mF. (7) 

ic 


The Cuboc order is basically detectable by the 12- 
sublattice magnetization Map- As was mentioned in the 
previous section, moreover, there is another essential or¬ 
der parameter —chirality degrees of freedom— originat¬ 
ing from the 120° structure on the unit triangle. For the 
unit triangle at r, we define the vector spin chirality^, 

«^(^) = ^ ['S'a(t') X Ssir) + Ssir) x 5'c(r) 

-i-Sc{r)x SA{r)]. ( 8 ) 

An important point on the Cuboc order is that the vec¬ 
tor chiralities associated with the four triangles in the 
extended unit cell also have non-coplanar configuration. 
Thus, we introduce the sublattice vector chirality order 
parameter for the Cuboc state as 

(9) 

^ ic ^13 

where J2rp ' sums up k, carrying the sublattice index /3 
of the extended unit cell. Here, note that the vector 
spin chirality takes the same sign in the c-axis direc¬ 
tion, although the direction of the spins are alternatingly 
aligned. 

In the Cuboc phase, both of Ma,p and Kp have finite 
values, but Mp = 0. In the ferromagnetic state, Mp 
is finite, while Kp = 0. For the incommensurate phase 
between the Cuboc and ferromagnetic phases, we do not 
set up a direct order parameter of the incommensurate 
spin configuration, because it is very difficult to deter¬ 
mine the pitch of the incommensurate oscillation within 
the system size up to L = 36. However, it should be 
noted that Ma,p = 0 and Kp = Q\s basically seen in the 
incommensurate phase. 


IV. RESULTS 

We present results of MC simulations for typical pa¬ 
rameters of the coupled spin tubes. As mentioned in Sec. 
H, the parameters corresponding to the CsCrF 4 is that 
J 1 /J 2 ~ 2 and J 3 /J 2 — 0. In this paper, we fix J 2 = 1.0 
and investigate the inter-tube-coupling (J 3 ) dependence 
for Ji = 1.0 and 3.0, along the arrow with the dotted line 
in the ground-state phase diagram as in Fig. [2{b). We 
then find no qualitative difference between Ji = 1.0 and 
3.0, so that we show the results for Ji = 1.0 below. 


A. Phase diagram 

In Fig. 01 we show the final phase diagram in the T-J 3 
plane, before presenting detailed analysis of simulations. 
Note that the horizontal axis represents — J 3 , since we 
consider the negative J 3 region. In the figure, it is verihed 
that the long-range orders at a finite temperature are 
consistent with the ground-state phase diagram of the 
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planar Kagome-triangular model in Fig. [2j For 0 < 
— J 3 < 1, the Cuboc phase is actually realized, and, in 
1 < — J 3 < 2 , the incommensurate spin order appears. 
Moreover, the system exhibits the ferromagnetic order in 
—J 3 > 2 , where the ferromagnetic coupling is dominant. 



FIG. 4: T-Jz phase diagram of the coupled triangular spin 
tubes with Ji = J 2 = 1.0. The Cuboc order is realized at 
finite temperature by the inter-tube interaction. 
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FIG. 5: (Color online) Results of MC simulations for 

— J 3 = 0.5. (a) Specific heat C. (b) Mean-square average 
of the Cuboc sublattice magnetization {Ma,i^)- (c) Binder 
cumulant for Ma,i- (d) Finite-size-scaling plot for the Cuboc 
order parameter, which yields Tc = 0.64721(4), u = 0.435(4) 
and 7 = 0.90(2). 


The phase boundary of the Cuboc-disorder transition 
is of second order for — J 3 < 0.7, which can be determined 
with a finite-size-scaling analysis. The universality of this 
second-order-transition line possibly belongs to a novel 
universality class associated with the triple-q structure of 
the Cuboc order. A detailed scaling analysis will be pre¬ 
sented in the following subsection. For 0.85 < — J 3 < 1.1, 
however, the Cuboc-disorder transition changes to the 
first order, where the double peak of the energy his¬ 
togram is observed. Note that a precise identification 
of the transition is difficult around — J 3 ~ 0.8, where a 
tricritical point is expected. In 1.1 < —J 3 < 1.8, the 
incommensurate-disorder transition line is estimated by 
a peak position of specific heat within finite size systems 
up to A = 32, where the size extrapolation is difficult. For 
—J 3 > 1 . 8 , the line of the disorder-ferromagnetic order 
transition is of second order, for which the universality 
class is consistent with the 3D ferromagnetic Heisenberg 
model. On the other hand, the transitions between the 
incommensurate phase and the other ordered phases are 
expected to be a commensurate-incommensurate type. 
These transition lines associated with the incommensu¬ 
rate order are also estimated within finite size results. 


B. Cuboc phase 

Let us begin with detailed analysis of the Cuboc phase. 
In Figs. [5]and[6l we present MC results for the Cuboc 
phase. Specific heat C and a mean-square average of 
the Cuboc sublattice magnetization (AFa,!^) are shown 
in Fig. [SKa) and (b). First of all, the specific heat C 
has a sharp single peak at Tc ~ 0.65, and the sublat¬ 
tice magnetization Ma.i also exhibits the phase tran¬ 


sition behavior at the same temperature. In order to 
precisely determine the transition point, we calculate the 
Binder cumulant^ of the Cuboc sublattice magnetization 
(Ma,i^)/(Ma,i^)^. The result for — J 3 = 0.5 is shown in 
Fig. [5](c), where the curves for various system sizes cross 
at Tc = 0.6470(5). This implies that the transition is of 
second order, prompting us to determine the universal¬ 
ity class of the Cuboc transition with a finite-size-scaling 
analysis. Assuming the scaling form for the susceptibility 
of the Cuboc sublattice magnetization, 

X = LNp{MA,i^)/T cx ( 10 ) 

with t = (T—Tc)/Tc, we perform the Bayesian estimation 
for the critical exponents ly, 7 and Tc^ Fig. [SKd) shows 
the resulting finite-size-scaling plot for y with the best- 
fit value V = 0.435(4), 7 = 0.90(2) and Tc = 0.64721(4). 
Here, we note that Tc is consistent with the result of the 
Binder cumulant, although it is obtained independently 
of the Binder-cumulant result. Taking account of the 
error originating from choice of the data window, we fi¬ 
nally adopt v = 0.44(2) and 7 = 0.91(3) for the critical 
exponents of the Cuboc sublattice magnetization. 

As mentioned before, the sublattice vector spin chi¬ 
rality is also another essential order parameter of the 
Cuboc order. Fig. [5] (a) shows {Ki^) for — J 3 = 0.5, 
where the transition occurs at the same Tc as the sub¬ 
lattice magnetization Ma.i- This behavior is consistent 
with the observation that the specific heat has a sin¬ 
gle peak at Tc. Thus, we can expect that the Cuboc 
magnetization and the chirality degrees of freedom ex¬ 
hibit the simultaneous transition. We then perform the 
finite-size-scaling analysis for the chirality susceptibil¬ 
ity with xk = TNp{Ki'^)/T oc In 
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FIG. 6 : (Color online) Results of the vector spin chirality for 
— Js = 0.5. (a) Mean-square average of the sublattice vector 
spin chirality (Ki^). (b) Finite-size-scaling plot for 
which yields best fit values: Tc — 0.64742(4), u = 0.433(8) 
and 7 x = 0.57(2). 


Fig. [S] (b), we present the finite-size-scaling plot with 
Tc = 0.64742(4), jz = 0.433(8) and -fK = 0.57(2), which 
are also obtained with the Bayesian estimation^. Taking 
account of the data-window dependence, we finally iden¬ 
tify the exponents associated with the vector spin chiral¬ 
ity BS V = 0.43(2) and jk = 0.55(4). Here, it should be 
noted that, although no apriori assumption of Tc and v 
was set up in this scaling analysis, the resulting Tc and 
v are consistent with those for y, while ■jk is clearly dif¬ 
ferent from 7. These facts suggest that the singular part 
of the free energy is scaled with 

fs cx hLy, hKLy^) ( 11 ) 

where y(= (^ -b d)/ 2 ) and ?/k(= + d)/2) are 

the eigenvalues of the linearized-renormalization-group 
transformation corresponding to the fields conjugated 
to the Cuboc magnetization and the chirality, respec¬ 
tively. Thus we have concluded that the transition of 
the spin and chirality degrees of freedom is simulta¬ 
neous. Here, the exponents obtained for the Cuboc 
transition are clearly different from those of the lay¬ 
ered triangular lattice antiferromagnet, although the 
simultaneous transition was also observed^^^— . The 
universality of the Cuboc transition provides a novel 
class, which might be characterized by the effective chi¬ 
ral Ginzburg-Landau-Wilson theory associated with the 
triple q-structure2^i21i2S. However, the detailed analysis 
of the effective model is an interesting future issue. Here, 
we remark that Mci,p and Kp are confirmed to be con¬ 
sistent among the all combinations of sublattice indices 
a and /3. On the basis of the analysis above, we have fi¬ 
nally drawn the second-order-transition line of the Cuboc 
phase in Fig. 01 

We turn to the hrst-order transition in the region of 
0.85 < — J 3 < 1.1. Figure [7] shows MC results for — J 3 = 
1.0. In Fig. [7Ka), the specific heat C also has a sharp 
single peak at Tc 0.675. In Fig. [T^b) and (c), moreover, 
we can observe that mean-square averages of the Cuboc 
sublattice magnetization (Ma,i^) and of the vector spin 
chirality {Ki ) also exhibit the phase transition behavior 
at the same temperature Tc, illustrating the simultaneous 
transition of the spin and chirality degrees of freedom. In 



FIG. 7: (Golor online) Results of MG simulations at — J 3 = 
1.0. (a) Specific heat C. (b) Mean-square average of the 

Guboc sublattice magnetization {Ma,i‘^)- (c) Mean-square 
average of the sublattice vector spin chirality (d) En¬ 

ergy histogram around the transition point for L = 32. The 
double peak structure emerges at Tc = 0.675, indicating that 
the transition is of first-order. 


the figures, we can also observe that the both of {Ma,i^) 
and {Ki^) for L = 24 and 32 show small jumps at Tc, 
suggesting that the transition is of first order. We have 
computed the energy histogram around Tc to confirm its 
double peak structure at Tc = 0.675 (Fig. E^d)). Thus, 
we conclude that the transition in 0.85 < — J 3 < 1.1 is 
the hrst-order. 

Finally, we would like to comment on the tricritical 
point expected around — J 3 ~ 0.8. As varying J 3 , we have 
checked that the double peak of the energy histogram 
appears down to —J 3 = 0.85, while the crossing point 
of the Binder cumulant emerges up to — J 3 = 0.7. Thus, 
the tricritical point is possibly located around — J 3 ^ 0 . 8 . 
Within the present system size of the MC simulation, 
however, a precise identihcation of the tricritical point is 
difficult. The detailed analysis on this respect is a future 
issue. 


C. Ferromagnetic phase 

In the — J 3 > 2.0 region, the ferromagnetic coupling 
becomes dominant, where the system forms the ferro¬ 
magnetic order in a Kagome-triangular layer, accompa¬ 
nying the second-order phase transition. In Fig. 01 we 
show results of MC simulations at — J 3 = 3.0. The spe¬ 
cific heat C in Fig. Ha) indicates a peak at Tc ~ 2 . 2 , 
and the magnetization (AFp^) in Fig. Hb) also exhibits 
the phase transition behavior. Note that the sublattice 
spin chirality {Ki^) is checked to be always zero in the 
region of the ferromagnetic phase. 

In order to precisely determine Tc, we further calcu¬ 
late the Binder cumulant {Mp'^)/. Its result is 
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FIG. 8 : (Color online) Results for — J 3 = 3.0. (a) Specific 
heat C. (b) Mean-square average of the uniform magnetiza¬ 
tion (AFf^). (c) Binder cumulant for My. (d) Finite-size¬ 
scaling plot for the susceptibility of xf- The fitting result is 
consistent with the 3D Heisenberg universality class^^. 


shown in Fig. He), where the crossing point appears at 
Tc = 2.233(3). We also perform the finite-size-scaling 
analysis of the susceptibility xf = ^‘^LNp{Mp‘^)/T, us¬ 
ing the Bayesian estimation. The finite-size-scaling plot 
in Fig. Hd) well collapses to a scaling function with 
Tc = 2.235(5), 1 / = 0.70(2) and 7 = 1.40(5). Note that 
these exponents are clearly consistent with the 3D ferro¬ 
magnetic Heisenberg class: ly ~ 0.71 and 7 ~ 1.40^“—. 
Thus, we can conclude that the universality class of the 
transition for the ferromagnetic order is the conventional 
3D ferromagnetic Heisenberg class. 


D. Incommensurate phase 

Let us hnally discuss the incommensurate phase in the 
1.1 < — J 3 < 2 region, which is sandwiched between the 
Cuboc and ferromagnetic phases. In Fig. we show 
the specific heat (7 at — J 3 = 1.5, where a phase transition 
is illustrated by a rounded peak of C around Tc ^ 0.7. 
However, the size extrapolation to extract the bulk be¬ 
havior is usually difficult for the incommensurate order, 
where the pitch of the oscillation does not match with the 
system size. Thus, we basically estimate the transition 
temperature by the peak position of the specific heat of 
L = 32. 

Turning to the transition line between the incom¬ 
mensurate and ferromagnetic phases, we can expect the 
commensurate-incommensurate type transition, where 
the wave vector of the order may continuously sift from 
the F point toward the M point, following the ground- 
state phase diagram in Fig. 2(b). Figure Hb) shows the 
uniform magnetization {Mp^) for — J 3 = 1.9, where we 
observe that it takes a finite value in 0.35 ^ T < 1.2, 
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FIG. 9: (Color online) Results of MG simulations for the 
incommensurate phase, (a) Specific heat C for —J 3 = 1.5, 
which shows a round peak at T ~ 0.7. (b) Mean-square aver¬ 
age of the ferromagnetic magnetization {Mp^) at — J 3 = 1.9, 
which abruptly decays in the incommensurate order region. 


but rapidly decays below T ^ 0.35. This behavior in¬ 
dicates that the ferromagnetic order appears between 
0.35 ^ T < 1.2, but it abruptly changes into the in¬ 
commensurate order in the low-temperature region of 
T < 0.35. We thus define the boundary between fer¬ 
romagnetic and incommensurate phases as the middle 
point of the onset and offset of (ATp^) for L = 32. 

On the other hand, we note that the transition line 
between the incommensurate and Cuboc phases is diffi¬ 
cult to estimate by the result within L = 32. Thus, the 
border between the incommensurate and Cuboc phases 
in Fig. Hlis just a guide for eyes. 


V. SUMMARY AND DISCUSSIONS 

In this paper, we have investigated phase transitions 
of the coupled triangular spin tubes associated with 
CsCrF 4 . An essential point is that a two-dimensional 
section of the coupled tubes forms a Kagome-triangular 
plane, which drives the system into exotic orders such as 
Cuboc order, incommensurate order and ferromagnetic 
order. In particular, the Cuboc order is characterized 
by the twelve sublattice non-coplanar spin structure car¬ 
rying the triple-q wave vectors, which accompanies the 
non-coplanar structure of the vector spin chirality as well. 
Performing extensive Monte Carlo simulations, we have 
demonstrated that these phases are actually realized for 
the coupled tubes with the negative inter-tube coupling 
J 3 at a finite temperature. The resulting phase diagram 
was summarized in Fig. IT] Then, a particular finding 
is that the transition to the Cuboc order in — J 3 < 0.7 
is described by the simultaneous second-order transition 
of the spin and chirality degrees of freedom. To our 
knowledge, the universality of this transition is a novel 
class characterized by the non-coplanar spin structure 
accompanying the chirality, which might be described 
by an effective chilal Ginzburg-Landau-Wilson theory 
associated with the triple g-structurei ^^7^i^^ While for 
0.85 < — J 3 < 1.1, the transition is of first order, where 
the double peak of the energy histogram is confirmed. 
However, the analysis of the expected tricritical point is 
a remaining issue. On the other hand, for — J 3 > 1.9, we 












have confirmed that the ferromagnetic transition belongs 
to the 3D Heisenberg universality class. 

From the experimental view point, CsCrF 4 is de¬ 
scribed by the weak J 3 coupling limit of the present 
model. The neutron scattering experiment of CsCrF 4 
actually suggests that a possible spin order is not a naive 
120 ° structure;^ and thus a finite temperature transi¬ 
tion to the Cuboc phase can be expected. However, 
we should also take account of another fact that the 
specific-heat experiment of CsCrF 4 captures no anomaly 
down to 1.5K, while the bulk phase transition to the 
Cuboc phase should theoretically accompany a certain 
anomaly of the specific heat. A reason for this incon¬ 
sistency is an anisotropic interaction effect. Since the 
inter-tube coupling J 3 of CsCrF 4 is basically very small, 
the Dzyaloshinsky-Moriya interactional^, which is actu¬ 
ally suggested in CsCrF 4 due to its crystal structure^^, 
can compete with the small J 3 coupling. Then, such an 
anisotropy effect may affect the Cuboc order configura¬ 
tion, and the transition with a small scale anomaly could 
be easily modified into a weak crossover. A direct com¬ 
parison of the neutron scattering experiment of CsCrF 4 
with the Cuboc order is highly desired. In addition to 
the above, we should also analyze how the quantum fluc¬ 
tuation affects the stability of the Cuboc order, which 
is another significant problem to understand the CsCrF 4 
experiment. 

In this paper, we have basically investigated the strong 
leg-coupling region (Ji > J 2 » 1 *^ 31 ), since our moti¬ 


vation is in the spin tube system. Our result implies 
that the leg-coupling certainly stabilizes the Cuboc or¬ 
der to be a true long-range order with the simultane¬ 
ous transition of the spin and chirality. On the other 
hand, the Cuboc order was originally proposed for the 
planar Kagome model with the next-nearest-neighbor 
coupling, where the spin and Zi chirality transitions are 
separated2ii2^. Recently, a Cuboc-type spin fluctuation 
was actually suggested for Kapellasite, which may be de¬ 
scribed by a S' = 1/2 Kagome antiferromagnet contain¬ 
ing up to third-nearest-neighbor couplings^. Theoreti¬ 
cal investigations stimulated by Kapellasite also have re¬ 
vealed interesting properties attributed to the Kagome 
structure4^“— Thus, it is an essential problem to under¬ 
stand how the 3D Cuboc class can be connected to the 
spin-liquid-like behavior with the Z 2 -chirality breaking 
in the limit of the planar Kagome-triangular model. 
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